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A  MATHEMATICAL  MODEL  FOR  AIR  FLOW  IN  A  VEGETATIVE  CANOPY 


RONALD  M.  CIONCO  I 
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INTRODUCTION  , 
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Ttio  usual  concept  of  _,urfaco  roup*hne*\s  Anri  h" 
5f~°fpJ?ne  displacement  breaks  dov/n  when  one  ccnsidors  a" 
vegetative  canopy  as  the  lower  boundary  for  the  planotarv 

?Q  ?ouShnoss  of  many^egetati  ve ^canopies 
fr°m  the  log  law  changes  with  wind  speed  as  do  the 
^SC««a?K  <UsP15oeme"ta>  Ideally  one  would  like  to  " 

ofPfte  ho!i‘L?0rSdyn?1?i0  r°uehness  of  vegetation  In  terras 
of  Its  height ,  density  and  drag  characteristics , 

Restricting  considerations  to  the  turbulent 
momentum>  a  model  was  developed  that  will 
predict  the  canopy  wind  profile  within  semi-rii/id  canonic 
For  the  purpose  of  this  report,  a  oanony  Is  So fined  as 
V.hat  layer  spanning  the  region  from  the  ground  surface  to 
the  top  of  the' plant.  To  date  three  canopy  flow  models 

third  moSs?eVeT?PS?'  ^iS  pap9r  Pri“rllyy describes  the 
third  model.  It  also  discusses  the  two  previous  models 

ofmanre'  ?^atsd°  and  Appleby,  1963)  briefly,  and  then 
compares  this  prior  work  to  that  of  other  instigators . ' 

2.  BACKGROUND 

The  earliest  efforts  in  canoov  modelling 

RlttPred(,bLLeni0n  iaSd  Sf?ller>  Tan  and*  Ling,  and  Ordway. 
Ritter,  upence  and  Tan  (Lemon  et  al,  1963),  In  lc6o  * 

r  way  et  al  (I960)  proposed  a  simple  canopy  flow  model 

1 1  was  suggested  that  the  turbulent  transfer  In  the  mL, 

could  be  pursued  In  a  manner  similar  to  that  of  the  "  Py 

.surface  boundary  layer,  provided  allowance  is  made  for  the 

loss  of  momentum  to  the  leaves  and  stalks  of  the  vegetal 

tion.  They  further  proposed  that  the  loss  of  momentum  wa- 

proportional  to  the  square  of  the  local  mean  velocity. 
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This  condition  is  analogous  to  the  definition  of  the  drag 
coefficient  for  fully  rough  flow  through  a  pipe.  Under 
this  condition  the  drag  coefficient  is  independent  of  the 
Reynolds  number.  Also,  for  the  steady  state  condition 
with  no  advection,  the  local  loss  of  momentum  must  be 
equal  to  the  convergence  of  momentum  transport.  To 
simplify  their  model,  Ordway  et  al  assumed  the  transfer 
coefficient,  Km*  to  be  constant  within  the  canopy.  In  a 
later  report,  Tan  and  Ling  (1961)  allowed  Km  to  increase 
linearly  with  height  as  it  does  in  the  surface  boundary 
layer.  The  author’s  initial  canopy  flow  model  (Cionco, 
1962)  is  similar  to  these  models  in  many  respects  with 
one  Important  exception j  the  model  does  not  assume  what 
the  behavior  of  Km  is  within  the  canopy,  and  rather  than 
working  with  the  transfer  coefficient,  it  was  more  desir¬ 
able  to  investigate  the  mixing  length  and  its  properties, 
Thus  for  turbulent  transfer  of  momentum  within  vegetative 
canopies  the  following  equations  result: 


or 


ix 

3  z 


pSu2 


(1) 
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(2) 


whera  t  is  the  shearing  stress,  z  is  the  vertical  space 
coordinate,  p  is  the  density  and  u  is  the  mean  wind 
velocity,  S  is  defined  as: 


S-icfcA  (3) 

and  is  similar  to  Prandtl’s  skin  friction  coefficient, 

CA  is  the  drag  coefficient  of  the  canopy  elements,  A  is 
the  effective  aerodynamic  surface  area  of  the  vegetation 
per  unit  volume  and  l  is  the  mixing  length.  This  model 
yielded  a  mixing  length  that  was  essentially  constant 
with  height  in  a  corn  plant  canopy.  The  wind  profile 
solutions  for  a  semi-rigid  canopy  using  these  general 
mixing  characteristics  in  equation  (2)  were  in  excellent 
agreement  with  the  observed  canopy  wind  profiles, 


The  second  model  formulated  a  canopy  of  elements 
that  are  ideally  uniform  in  geometry,  distribution  and 
behavior.  The  basic  assumptions  of  this  "ideal"  canopy 
model  are  that  the  scale  length  of  turbulence  and  the 
dissipation  coefficient  are  constant  with  height.  Thus 
the  idealised  canopy  can  be  defined  as  having  (1)  a 
uniform  vertical  leaf  area  distribution,  (2)  a  uniform 
vertical  distribution  of  the  drag  coefficient  of  the 
canopy  elements,  and  (3)  a  drag  coefficient  that  is 
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independent  of  the  local  Reynolds  number*  The  properties 
of  this  model  are  characterized  by  (1)  a  mixing  length 
that  is  constant  with  height  except  very  near  the  surface, 

(2)  the  turbulence  intensity  is  constant  with  height,  and 

(3)  the  wind  velocity  is  an  exponential  function  of  height 
of  the  following  form: 

ft  a(J  -  1) 

U  «  UH  e  H  ( in 


where  the  attenuation  coefficient  "a”  is  a  constant.  The 
model  also  defines  relationships  of  the  canopy’s  basic 
parameters  in  terms  of  "a”  in  the  following  forms: 


and 


(5) 

(6) 


where  ”a”  is  constant  if  both  S  and  a0  are  unchanging 
with  height.  A  third  relation  evolves  from  the  above 
a-terms  as: 


m 

2L20 


(7) 


It  relates  the  drag  at  the  oanopy  top  to  the  mixing 
length  within  the  canopy. 


When  applying  the  ’’ideal”  model  to  the  natural 
canopy  during  the  steady  state  condition  with  no  advoc- 
tion,  it  was  found  that  (1)  the  wind  profiles  within  a 
corn  plant  oanopy  reported  by  Tan  and  Ling  (1961)  were 
exponential  and  that  the  mixing  length  was  computed  to 
be  constant  with  height,  (2)  results  reported  by  Makagawa 
(1956)  for  a  rice  crop  and  by  Inoue,  Lemon  and  Dcnmoad 
(196*0  for  a  wheat  field  have  shown  turbulence  Intensity 
to  be  almost  independent  of  height,  and  (3)  limber  crops 
like  alfalfa  (Lemon  et  al,  1963)  wore  definitely  not 
*  ideal  in  nature.  Cionco  et  al  also  found  that  S,  Cp,  and 
Ac  were  not  conservative  properties  of  the,  canopy,  but 
were  interdependent  and  only  ”a"  as  expressed  in  equation 

(5)  was  nearly  constant.  The  other  relations  in  equations 

(6) '  and  (7)  were  sensitive  to  windspeed,  Salto  (196*0 
also  points  out  that  the  condition  ”a”  «  constant  through 
all  heights  does  not  seem  to  be  satisfied  in  the  general 
case « 
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Other  investigators  working  with  more  recent 
canopy  data  have  found  similar  results.  Inoue  (1963) 
discussed  corn,  rice  and  wheat  canopies.  Inoue,  Lemon, 
and  Denmead  (196*1)  reported  a  wheat  canopy,  Uchijima  and 
Wright  (1963)  described  an  immature  corn  canopy,  and 
Salto  (196**)  discussed  a  wheat  canopy  and  two  separate 
immature  corn  canopies. 

Inouo,  Saito,  and  Inoue  ot  al  have  also  found 
that  the  constant  mixing  length  and  the  exponential  wind 
equation  apply  to  their  canopy  studies.  Except  for  some 
of  the  immature  crops,  their  attentuation  coefficients 
were  essentially  const  ant  for  the  respective  crop  canopies. 
Uchijima  and  Wright  fitted  a  log-linear  function  of  height 
and  velocity  for  their  profiles;  they  also  described 
the  mixing  length  to  increase  linearly  with  height. 

Howevor,  Inoue  (1963)  pointed  out  that  the  data  does  not 
support  their  conclusion  for  i.  He  analyzed  1  to  bo 
essentially  constant  throughout  most  of  the  canopy.  One 
other  interesting  result  of  these  analyses  is  that  Inoue 
ot  al  and  Inoue  also  established  the  relation  proposed  by 
Cionco  et  al  shown  as  equation  (6). 

3.  MIXING  LENGTH  MODEL  AND  SOLUTIONS 

The  studies  and  models  described  thus  far  all 
define  the  resistive  effect  of  the  vegotation  Jn  the  form 
proposed  by  Ordway  et  al  (I960).  They  follow  Prandtl*s 
skin  friction  concept  (Prandtl  and  Tietjens,  1957). 

Howqver,  the  new  model  will  depart  from  this  concept  to 
ptfesont  a  more  realistic  simulation  of  the  natural  canopy. 
This  change  ia  suggested  by  the  work  reported  by  Kutabaoh 
(1961),  Empirical  results  derived  from  his  field  studies 
indicated  anothor  form  of  S  (defined  by  equation  (3)) 
would  be  more  appropriate  for  field  observations. 

Kutzbach  conducted  an  experiment  to  determine  the  effect 
of  artificial  roughness  elements  upon  the  wind  profile 
over  an  ice  covered  lake.  The  roughness  elements  were 
ordinary  bushel  baskets  and  the  experiment  investigated 
the  orderly  variations  of  both  the  density  and  heights 
of  the  roughness  elements.  The  data  revealed  that  3 
should  be  redefined  so  that  the  dens ity ' function  of  the 
roughness  elements  has  more  meaning.  Further  analysis 
of  his  observations  indicated  that  the  parameter  A  should 
be  a  squared  term,1  To  keep  S  dimensionally  correct  , and 

‘The  bushel  basket  analysis  yielded  C^aCp0*7  and  also 
CpaA1*243  so  that  SaA1*0,  For  computational  convenience, 
1,8  was  rounded  off  to  2  so  that  3aA2,  See  page  25  of 
Cionco  et  al  (1963)  for  these  relationships , 
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consistent  with  the  previous  theoretical  model  (Cionco  ot 
al,  1963)  an  additional  length  term  is  necessary.  The 
most  obvious  length  parameter  to  be  considered  for  the 
canopy  flow  model  is  the  mixing  length  in  the  crop. 
Equation  (8)  is  the  definition  of  S’  based  on  the  bushel 
basket  results. 


S'  «  A  0Jk  A*t  (8) 

2  D 

Equation  (2)  can  now  be  written  making  use  of 
equation  ( 8) : 

h[lT%Y  ») 

A  solution  of  the  above  equation,  of  course,  poses  a 
problem  in  that  there  are  more  unknowns  than  equations. 

All  of  the  components  of  S’  are  not  known.  Cft  in  partic¬ 
ular  is  unknown.  A  is  known  and  can  bo  evaluated 
graphically  from  the  agronomist’s  cumulative  leaf  area 
index  shown  in  Figure  1  (Allen  et  al,  196*1).  The  mixing 
length  in  also  unknown  and,  of  course,  the  velocity 
distribution  is  the  intonded  solution,  Of  these  unknowns, 
the  only  data  available  are  wind  profilos  in  the  canopy. 
This  fortunately  is  sufficient  information  to  allow  for 
fi  solution  of  ft  and  CJ  as  one  unknown.  The  intention 
now  is  to  provide  enough  information  to  establish  some 
general  properties  of  the  canopy’s  mixing  processes  and 
drag  characteristics. 


The  mixing-length  solution  of  equation  (9)  oan 
be  expressed  at  follows  in  a  normalized' form! 


L.  t  Jl.  1  m  a?y_2  -  f  h  1  a2Vax2 

ax  [  QLi/2j  sc  avr/ax) 2  [01,1/2]  T77T^  do) 


For  computational  purposes  It  was  desirable  to  express 
the  solution  in  nondimonsional  terms.  The  parameters  were 
arbitrarily  normalized  with  respect  to  a  level  twice  the 
height  of  the  canopy.  The  nondimonsional  variables  are 
X  *  Z/2H,  Y  ■  u/u2ji  and  L  »  ft2/fc2M,  The  D-tcrm  is  a 
collection  of  constants  along  with  Cl  and  is  defined  as 
S  ■  8H3(*  cn)/y-2u»  Equation  (10)  is  simply  the  result  of 
differentiating  equation  (9)  by  parts  arid  dividing  through 
by  L't  Thus  given  a  canopy  wind  profile  and  the  vertical 
leaf  area  distribution,  equation  (10)  will  yield  a  solu¬ 
tion  of  L  and  CJ, 
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The  computer  form  of  equation  (10)  is: 

H  “  C‘  +  EC*  (11) 

where  E  *  L  B"1  L~l//2  and  Cl  and  C2  are  the  respective 
righthand  terms  of  the  previous  expression.  Equation  (11) 
is  concerned  with  two  distinct  regions  of  turbulent  flow— 
that  within  the  canopy  and  that  immediately  above.  For 
Z>H|  Ci  must  be  zero  for  this  expression  to  apply  in  the 
region  where  the  logarithmic  wind  profile  is  valid, 

The  first  mixing  length  solutions  of  equation 
(11),  wero  calculated  utilizing  Tan  and  Ling*s  (1961)  two 
complete  wind  profiles  shown  in  Figure  2,  These  same 
complete  profiles  are  shown  in  Figure  3  to  depict  their 
fit  of  the  exponential  wind  profile  of  equation  ( *0  , 

Second  and  third  or<ler  lagrangian  polynomials  were  used 
to  estimate  tho  velocity  derivatives  of  Ci  and  C2  from 
the  normalized  data,  Tho  initial  conditions  for  the 
,2  "  0  level  are  Cj  ■  zero  and  C,  defined  by  2H/zoa  where 
z0  was  tho  roughness  length  at  the  ground  surface,  E  os 
defined  for  equation  (13)  was  approximately  10- 4  for 
2*0  and  B  is  of  the  order  of  102  for  a  semi-rigid 
canopy  crop,  The  A-torm  was  constant  from  the  ground 
surface  to  tho  uppermost  2055  in  which  tho  donsity  of  the 
foliage  docroaood  uniformly  to  zero  at  the  top  of  the 
vegetation,  Equation  (11)  was  then  integrated  from  the 
ground  surface  to  2H  using  a  fourth-order  Runge-Kutta 
integration  routine  ,(Niolfiori#  1956), 

The  solutions  are  shown  in  Figure  4  as  profiles 
of,  the  normalized  mixing-length  term,  E,  versus  the 
normalized  height,  X,  Tho  height  has  been  normalized  so 
that  tho  velocities  and  mixing  lengths  would  bu  repre¬ 
sentative  of  a  unit  canopy.  Within  tho  canopy,  the  mixing 
length  is  essentially  constant  except  near  the  ground 
surface i  above  the  canopy  the  mixing  length  increases 
linearly  with  height.  The  solution  for  Corn  1  firmly 
supports  these  conclusions,  while  the  results  for  Corn  2 
are  less  consistent  within  the  canopy,  • 

In  a  general  way,  these  solutions  are  compatible 
with  the  idealized  canopy  concept  of  Cionco  et  al  and 
conform  to  the  canopy  scheme  proposed  by  Inoue,  Tho 
nearly  constant  mixing  lengths  calculated  from  the  pro¬ 
files  confirm  the  interdependence  of  S’  and  l  with  height, 
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4,  CANOPY  FLOW  MODEL  AND  SOLUTIONS 


In  the  previous  section,  it  was  necessary  to 
establish  some  general  characteristics  of  the  canopy’s 
mixing  processes.  With  this  information  now  available, 
a  solution  of  the  canopy  wind  profile  is  possible. 
Equation  (9)  can  be  expressed  as: 
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au  f  du 
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for  Z<H,  All  parameters  except  X  retain  the  same  defini¬ 
tions*^  previously.  The  model  treats  the  canopy  as  a 
unit,  thus  X  ■  Z/H, 


The  mixing  length  solutions  also  confirm  Inouo’s 
(1963)  representation  of  the  canopy  mixing  length  by  two 
expressions:  ft  ■  kz  near  the  earth’s  surface  and  ft  • 
constant  throughout  the  rest  of  the  canopy,  By  aooooting 
this  generalized  mixing  length  in  the  canopy,  3(lnft0)/3X 
oan  bo  determined.  To  allow  for  the  above  changes  in  the 
mixing  length  in  the  righthand  term,  a  generalized  ft2/ft$ 
as  a  function  of  hoight  was  oomputod  from  the  mixing 
length  solutions.  Defining  two  new  variables  of  the 
righthand  term  as  C  *  H3S’/ft|  and  0  “ft2/ftg A2,  equation  (12) 
oan  now  be  written  as: 
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The  C-term  is  now  the  only  unknown,  'To  deter¬ 
mine  0,  an  empirical  relationship  is  necessary,  From  the 
ideal  canopy  concept,  several  relationships  were  establi¬ 
shed  as  equations  (5,  6,  and  7)  that  should  be  unique 
functions  of  the  canopy.  Equation  (5)  evaluated  from  tho 
mature  corn  data  proved  to  be  reasonably  constant  and 
equal  to  2,2,  Equation  (7)»  however,  varied  as  a  function 
of  wind  speed,  When  modiflod  to  0  ■  U^  H2  Cn/2 ft|  it 
became  reasonably  constant  and  oqual  to  713.5*  From  these 
relationships ,  C  can  easily  be  determined  as  a  function 
of  where  C  *  S&B/Uji* 


Tho  canopy  wind  profile  solution  is  obtained 
from  a  nonlinear  overrelaxation  algorithm  described  by 
Lieberstoin  (1959).  It  is  an  iteration  method  that 
continually  solves  for  the  best  estimate  (BIST)  of  the 
velocity  such  that  equation  (13)  converges  to  zero}  i.e., 
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Second  order  lagrangian  polynomials  are  used  to  evaluate 
the  first  and  second  derivatives  of  the  velocity.  The 
initial  conditions  of  the  algorithm  are  the  velocity,  U0| 
at  the  ground  surface,  the  velocity,  U^,  at  the  canopv/alr 
interface,  and  a  single  value  of  C.  U0  is  of  court*© 
zero*  Uh  is  observed  or  computed  from  the  log  law  rela¬ 
tionships,  and  C  is  evaluated  directly  from  the  mixing 
length  solutions  or  from  the  a  and  8  relationships ,  From 
this,  one  sees  that  a  solution  is  easily  attainable  with 
the  minimum  information  of  Uh • 

Canopy  flow  solutions  of  the  original  Tan  and 
Ling  data  (Figure  2)  are  shown  in  Figure  5,  This  plot 
shows  good  agreement  between  the  observed  and  simulated 
profiles  of  Corn  HI  and  02*  The  other  simulated  profiles 
are  in  good  agreement  with  the  partially  observed  profiles 
of  the  upper  third  of  the  canopy.  The  largest  differences 
between  all  of  the  observed  and  simulated  values  are  555 
of  that  volooity  range, 

t 

5,  OTHER  RESULTS  AMD  DISCUSSION 


More  solutions  wore  obtained  as  an  independent 
chook  of  the  canopy  wind  profile  model,  Figures  6  and  7 
summarize  the  model’s  simulations  of  the  observations 
reported  by  Plate  and  Cerrnuk  <  1963)  and  Wright  (1963)* 
Figure  6  consists  of  two  separate  sots  of  data  from 
mature  corn  crop  canopies,  Five  profiles  are  for  Wright’s 
data;  the  other  datu  points  represent  five  profiles  for 
Tan  and  Ling’s  data,  Figure  7  compares  the  simulations 
to  Plate  and  Cermak’s  wind  tunnel  observations,  Those 
profiles  wore  collected  within  and  above  a  In  oh 
artificial  canopy  of  semi-rigid  plastic  strip  elements, 
This  plot  compares  ten  previously  simulated  solutions 
(for  mature  corn)  to  ten  wind  tunnel  profiles,  Generally 
speaking,  tho  computations  were  in  good  agreement  with  the 
observations , 

Statistically,  tho  agreement  can  bo  summarized 
as  follows.  Error  analysis  of  the  corn  canopy  solutions 
yioidod  an  average  absolute  error  of  11,7  cm/aoc  and  a 
•  root-mean-square-orror  of  1*1,8  cm/soc,  The  artificial 
canopy  solutions  have  an  average  absolute  error  of  12,7 
cm/sec  and  a  root-mean-square-error  of  15.1  cm/sec.  From 
this  analysis  and  Figure  7,  it  is  apparent  that  velocity 
predictions  of  100  cm/soc  or  less  are  subject  to  sizeable 
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error  for  the  artificial  canopy.  Fortunately,  this 
velocity  range  is  usually  restricted  to  the  lower  third 
of  the  canopy, 

The  under  prediction  of  the  low  velocity  range 
in  the  artificial  canopy  may  be  duo  to  the  bending  of  the 
plastic  strips.  As  mentioned  previously,  Saito  concluded 
that  although  “a"  appeared  constant  for  the  corn  canopy, 
it  would  not  be  true  in  the  general  case.  Because  the 
computed  profiles  wore  simulated  for  mature  corn,  "a"  was 
considered  constant  for  the  plastic  strips  rather  than  a 
function  of  velocity  and  flexibility,  This  would  also 
contribute  to  an  under  prediction  of  velocities  deep  in 
the  canopy. 

The  above  suggests  that  the  attenuation  coeffi¬ 
cient  may  be  a  useful  method  to  study  various  oanopies 
because  the  exponential  profile  has  boon  observed  in  quite 
a  variety  of  vegetations.  In  particular  the  effeots  of 
different  leaf  typos,  densities  and  degrees  of  rigidity 
should  bo  investigated.  This  coefficient  should  do 
rolatod  to  the  overall  canopy  density  and  flexibility, 
Sparse  vogotation  whioh  can  bo  readily  penetrated  by  the 
eddies  from  above  would  tend  to  have  small  attenuation 
coefficients,  These  coefficients  should  increase  as  the 
oanopy  donsity  inoreasos.  With  changes  in  flexibility 
of  the  vegetation,  one  would  expect  the  coefficient  to 
inoroaso  with  wind  speed  as  the  plant  streamlines  to  the 
wind,  More  intensive  investigation  of  "a"  from  wind 
profiles  in  canopies  of  known  densities  and  elasticities 
may  help  to  establish  the  aerodynamic  characteristics  of 
vegetation  by  simple  relationships, 

4y  supplementing  the  model  with  wind  profiles 
above  the  canopy  and  canopy  temperature  and  humidity 
observations,  the  model  should  prove  useful  in  studying 
the  canopy’s  energy  transport  and  diffusion  processes, 

In  describing  tho  sink  and  source  distributions  of  the 
energy  balance  components,  values  of  these  canopy  fluxes 
could  bo  approximated,  Tho  diffusion  of  matter  into  the 
canopy  could  also  bo  estimated  using  the  functional 
relationship  of  eddy  viscosity  to  the  mixing  length,  In 
this  light  the  generalised  mixing  length  discussed  pre¬ 
viously  will  take  on  more  importance, 

6,  SUMMARY 

Although  the  model  has  limitations  in  that  it 
has  been  applied  only  to  semi-rigid  vegetation,  it  does 
yield  reasonable  results  for  canopies  that  meet  the 
model’s  requirements.  The  computed  mixing  length 
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solutions  showed  that  l  was  nearly  constant  with  height 
in  the  canopy.  This  agreed  closely  to  the  l  properties 
proposed  previously  by  the  idealized  canopy  concept  and 
later  by  other  investigators.  These  solutions  also  showed 
that  l  increased  linearly  with  height  above  the  canopy  as 
it  should  in  the  logarithmic  wind  profile  domain.  Finally, 
the  computed  canopy  flow  solutions  were  in  good  agreement 
with  the  observations.  These  solutions  further  verified 
the  concept  that  the  canopy’s  mixing  length  profile  is 
essentially  constant. 
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Figure  3*  The  Logar¬ 
ithm  of  the  Wind  Speed 
as  a  Function  of  the 
Relative  Height  Within 
a  Mature  Corn  Canopy. 
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Figure  k ♦  The  Normalized  Mixing  Length 
Terms  as  a  Function  of  Normalized  Height 
Within  and  Above  a  Mature  Corn  Canopy 
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Figure  5.  A  Comparison  of  the  Observed 
Wind  Velocities  both  Within  and  Above  the 
Canopy  to  the  Calculations  Within  a  Mature 

Corn  Canopy 
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Figure  6.  A  Simple  Scatter  Diagram  Comparing 
Culouluted  and  Observed  Velocities  for  a 
Mature  Corn  Canopy.  The  curve  is  the  line  of 
perfect  prediction  and  does  not  represent  the 
data. 


Figure  7.  A  Simple  Scatter  Diagram  Comparing 
Calculated  and  Observed  Velocities  for  an 
Artificial  Canopy.  The  curve  is  the  line  of 
ge^fect  prediction  and  does  not  represent  the 
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